inboth<- read.csv("in_panel.csv")

#VPP (non interacted) is year FE as it's a 2-year DiD 
aqqq1<- feols(part~iswhiteperson*voucher_yr*incss2005|as.factor(id) , data=inboth)
aqqq3<- feols(part~republican*voucher_yr*incss2005|as.factor(id) , data=inboth)
aqqq6<- feols(part~income_perc*voucher_yr*incss2005|as.factor(id) , data=inboth)
library(modelsummary)
modelsummary(list(aqqq1, aqqq3, aqqq6), output = "latex")

